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We present a measurement of the differential cross section for tt events produced in pp collisions at 
\fs = 1.96 TeV as a function of the transverse momentum {pr) of the top quark. The selected events 
contain a high-pT lepton {£) , a large imbalance in pr , four or more jets with at least one candidate for 
a 6 jet, and correspond to 1 fb~^ of integrated luminosity recorded with the DO detector. Objects in 
the event are associated through a constrained kinematic fit to the tt ^ WbWb — >■ Ivbqq'b process. 
Results from next and next-to-next-to-leading-order perturbative QCD calculations agree with the 
measured differential cross section. Comparisons are also provided to predictions from Monte Carlo 
event generators using QCD calculations at different levels of precision. 

PACS numbers: 14.65.Ha,12.38.Qk,13.85.Qk 



The transverse momentum {pt) of top quarks in ti 
events provides a unique window on heavy-quark pro- 
duction at large momentum scales. In the standard 
model (SM), the lifetime of the top quark is far shorter 
than the characteristic hadron-formation time of quan- 
tum chromodynamics (QCD), which provides access to 



the properties and kinematics of a "bare" quark, such as 
mass, charge, spin, and p^, that are almost unaffected 
by bound-state formation or final-state interactions 
The top quark is unique in that it has a mass close to the 
scale of clectroweak symmetry breaking. Detailed studies 
of the properties of this bare quark beyond the measure- 
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ment of its total production rate, such as the measure- 
ment of its quantum numbers and of its couphngs to other 
SM particles, may indicate whether the top quark plays 
a privileged role in the symmetry breaking. Focusing on 
details of the tt production, measurements of differential 
cross sections in the tt system test perturbative QCD 
(pQCD) for heavy-quark production and can constrain 
potential new physics beyond the SM e.g., by mea- 
suring the transverse momentum of the top quark • 

In this Letter, we present a new measurement of the 
inclusive differential cross section for pp ^ tt + X pro- 
duction at -y/s = 1.96 TeV as a function of the pr of 
the top quark. The measurement is corrected for detec- 
tor efficiency, acceptance and resolution effects, making 
it possible to perform direct comparisons with different 
theoretical predictions. The data were acquired with the 
DO detector at the Fermilab Tevatron Collider and cor- 
respond to an integrated luminosity of ~ 1 fb^^. This 
measurement was performed in the i+jets decay channel 
of tt WbWb Iv + hh+ > 2 jets, where i repre- 
sents an e or /i from the decay of the W boson or from 
— >■ r — 7^ ^. The dependence of the cross section on 
the pt of the top quark was examined previously using 
« 100 pb"^ of Tevatron Run I data &i ^/s = 1.8 TeV 
where no deviations from the SM were reported. 

The DO detector [l2| is equipped with a 2 T solenoidal 
magnet surrounding silicon-microstrip and scintillating- 
fiber trackers. These are followed by electromagnetic 
(EM) and hadronic uranium/liquid argon calorimeters, 
and a muon spectrometer consisting of 1.8 T iron toroidal 
magnets and wire chambers and scintillation counters. 
Electrons are identified as track-matched energy clus- 
ters in the EM calorimeter. Muons are identified by 
matching tracks in the inner tracking detector with those 
in the muon spectrometer. Jets are reconstructed from 
calorimeter energies using the Run II iterative seed-based 
midpoint cone algorithm with a radius of 0.5 13|. Jets 
are identified as originating from a h quark using an artifi- 
cial neural network (6 NN) which combines several track- 
ing variables [l^ . Large missing transverse energy, 
(the negative of the vector sum of transverse energies 
of calorimeter cells, corrected for reconstructed muons) 
signifies the presence of an energetic neutrino. Events 
are selected using a three-level trigger system, which has 
access to tracking, calorimeter, and muon information, 
and assures that only events with the desired topology 
or with objects above certain energy thresholds are kept 
for further analysis. 

The analysis uses similar data samples, event selec- 
tion, and corrections as used in the inclusive ti — >■ £-|-jets 
cross-section measurements detailed in Ref. [lB|- Events 
accepted by lepton-|-jets triggers are subject to additional 
selection criteria including exactly one isolated lepton 
with PT > 20 GeV/c and > 4 jets with pr > 20 GeV/c 
and |?7| < 2.5 [16i]; at least one jet must have pt > 40 
GeV/c. At least one jet is also required to be tagged 



by the b NN algorithm. Additionally, we require > 
20 GeV (25 GeV) for the e-|-jets (^-|-jets) channel and 
electrons (muons) with \t]\ < 1.1 (2.0). 

Our measurement uses the ALPGEN [l^ event gener- 
ator, with PYTHIA [3] for parton showering, hadroniza- 
tion, and modeling of the underlying event, to simulate 
the inclusive tt signal. A pythia sample serves as a 
cross check. The CTEQ6L1 set of parton distribution 
functions (PDFs) [3] was used with a common factor- 
ization and renormalization scale set to /i = mt -I- p!^*"^ 
for mt = 170GeV/c^. Backgrounds are modeled with 
ALPGEN+PYTHIA for VF-l-jets and Z-|-jets production, 
PYTHIA for diboson (WW, WZ, and ZZ) production, 
and COMPHEP [2^ for single top-quark production. The 
detector response is simulated using GEANT 21|. The 
simulated tt signal is normalized to the cross section 
measured by a dedicated likelihood fit in the same fi- 
nal state using the same event selections (including the 
^-tagging requirement) and data as Ref. [15|, namely to 
8.46to;97 pb at a top-quark mass mt = 170GeV/c^ (in 
good agreement with the value extracted in this study 
by integrating the differential cross section). The dibo- 
son and single top-quark backgrounds are normalized to 
their SM predictions, .Z-l-jets to the prediction from next- 
to-leading-order (NLO) pQCD, and M^-l-jets such that 
the predicted number of events matches the data before 
applying b tagging. 

The small multijet background, in which a jet is 
misidentified as an isolated lepton, is non-negligible only 
in the e-fjets channel. Its rate is estimated from data 
using the large difference in the probability of electro- 
magnetic showers of real electrons or misidentified jets 
to satisfy the electron selection criteria. The details of 
the sample composition and the observed yields before 
and after requiring the jets to be tagged as 6-jet are pre- 
sented in Table H] 

The selection yields 145 and 141 events in the e-|-jets 
and /i-fjets decay channels, respectively. The measured 
ti signal fraction is 0.79, indicating that this sample is 
suitable for detailed studies of tt production. A con- 
strained kinematic fit to the tt final state, which takes 



TABLE I: Expected yields for signal and backgrounds samples 
and observed event counts in e-fjets and /i-|-jets channels. 



sample 


e-f-jets /i-fjets 


tt 




131 108 


W+jets 




10 15 


Z+iets 




3.0 3.1 


single top 
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summed prediction 
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ground uncertainty 
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data 
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FIG. 1: The reconstructed top-quark mass compared with 
expectation. Hashed areas represent statistical and jet energy 
calibration uncertainties on the prediction. 



into account the unreconstructed neutrino and finite ex- 
perimental resolution, is used to associate leptons and 



jets with individual top quarks |22,|23[. The fit assumes 
equal masses for the two reconstructed top quarks and 
the two reconstructed W boson masses are constrained 
to 80.4 GeV/c^. All possible permutations of objects 
needed to produce the it system are considered, and the 
solution of fitted leptonic and hadronic top-quark four- 
momenta with the smallest (the goodness of the fit) 
is selected for further analysis. The 6-jet assignment in- 
formation is not used in the selection of the best permu- 
tation to avoid the associated efficiency loss. The effects 
of possibly selecting a wrong permutation when choosing 
the one with the best are taken into account in the 
corrections of the measurement to the parton level. The 
solution with the best (second best) corresponds to 
the correct assignment of the quarks from the decay of 
the ti pair in 48% (17%) of events. 

The reconstructed top-quark mass (jrit) from the best 
fit in data, simulated tt signal, and background is shown 
in Fig. [TJ There is good agreement between the data and 
the sum of signal and background expectations in terms 
of the shape, resolution, and mean of the distribution in 
mt (x^/NDF = 1.28). The spectrum of the top quark 
(for leptonic and hadronic entries) in data, together with 
predicted signal and background, is shown in Fig. [2] for 
the best solution but now refitted with a top-quark mass 
fixed to 170 GeV/c^ (the value used in the inclusive cross 
section measurement [is']) to improve resolution. To ob- 
tain a background-subtracted data spectrum, the signal 



purity is fitted using signal and background contributions 
as a function of pt , and applied as a smooth multiplica- 
tive factor to the data. The result is the background- 
corrected distribution shown as a solid line in Fig. [3l 

The reconstructed "Pt spectrum is subsequently cor- 
rected for effects of finite experimental resolution, based 
on a regularized unfolding method 24, |25| using a mi- 



gration matrix between the reconstructed and parton -pT 
derived from simulation. The size of the pt bins was cho- 
sen based on the requirement that the purity (the fraction 
of parton-level events which are reconstructed in the cor- 
rect pt range) is > 50%, as shown in Table |lTl This also 
results in -pT bins which are larger than the experimental 
resolution for the top quark p-p. The correlation between 
reconstructed and correct pT is > 80%. Figure [3] com- 
pares the reconstructed and corrected results as a func- 
tion of the pt of the top quark. The dependence of the 
unfolding on the parton spectrum shape in the migration 
matrix is tested by reweighting the distribution with arbi- 
trary functions. Shape variations of « 20% induce 2 — 6% 
changes in the differential cross section. A correction for 
acceptance from the dependence of the spectrum on kine- 
matic restrictions of reconstructed quantities is applied 
to the unfolded distributions. 

The measured differential cross section as a function 
of the pt of the top quark (using for each event the two 
measurements obtained from the leptonic and hadronic 
top quark decays), dcr/dpT, is shown in Fig. |4] and tab- 
ulated in Table IHIl together with the NLO pQCD pre- 
diction 0, [i^l. The statistical uncertainties are esti- 
mated by performing 1000 pseudo-experiments where, 
in each experiment, the background-corrected spectrum 
is allowed to vary according to Poisson statistics and 
is then unfolded using the regularized migration matrix 
(Table nil . The largest experimental uncertainties affect- 
ing the shape of the pT distribution include jet energy 
calibration in data and in simulation (1.5 — 5.0%), jet 
reconstruction efficiency (0.7 — 3.5%), and jet energy res- 
olution (« 0.5%). The residual dependence of the un- 



TABLE II: The migration matrix between the reconstructed 
(rows) and parton (columns) top-quark pt derived from 
ALPGEN ti events passed through full detector simulation. The 
matrix indicates the fraction of events migrated from a given 
parton bin to the reconstructed bins. The binning used for 
correlating reconstructed and parton levels of are given at 
the left and top, respectively. Results in bold print are for 
diagonal terms. 



VT (GeV/c) 


0-45 


45-90 


90-140 


140-200 


200-300 


300-400 


0-45 


0.530 


0.162 


0.062 


0.020 


0.003 


0.000 


45-90 


0.344 


0.578 


0.227 


0.072 


0.021 


0.000 


90-140 


0.103 


0.228 


0.560 


0.223 


0.055 


0.031 


140-200 


0.019 


0.029 


0.145 


0.581 


0.232 


0.071 


200-300 


0.002 


0.002 


0.006 


0.103 


0.650 


0.363 


300-400 


0.000 


0.000 


0.000 


0.001 


0.038 


0.535 



6 



O 
> 

2,10 

■a 



■a 1 



10-^ 



10' 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



D0 



• data, 1 fb"^ 



tt 



I I Background 



50 100 150 200 250 300 350 400 



top quark p^[GeV/c] 



u 
> 

2,10 

■a 



■a 1 



10 



10- 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I. 

D0 data, 1 fb"^ 

— reconstructed 
• corrected 



■1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



50 100 150 200 250 300 350 400 



top quark p^ [GeV/c] 



FIG. 2: The pr spectrum of top quarks (two entries per event) 
compared with expectation. Hashed areas represent statisti- 
cal and jet energy calibration uncertainties on the prediction. 



folded result on the top-quark mass is 2 — 6% for mt 
in the 170-175 GeV/c^ range. This additional uncer- 
tainty does not need to be considered for comparisons 
with models in which is set to 170 GeV/c^. For the 
main background sources, W^/Z-|-jets, we have also con- 
sidered the variations of the background shape caused 
by uncertainties in the k-factors and in additional scale 
factors for heavy-flavour jets. Other systematic uncer- 
tainties account for uncertainties in the modeling of 
the signal, estimated from the difference between ALPGEN 
and PYTHIA, for uncertainties in the PDFs and in the b- 
quark fragmentation. The uncertainty on the integrated 
luminosity is 6.1%. The systematic uncertainties quoted 
in the following combine the uncertainty on the normal- 
ization (independent of px) with the shape-dependent 
systematics. The total correlated systematic uncertainty 
is 9.6% (including the uncertainty on luminosity) and 
the total systematic uncertainty on the cross section, in- 
tegrating over pt, is 10.7%. 

Results from NLO pQCD [2£ 
using CTEQ61 ^ PDFs (using the scale fi 
170GeV/c^) are overlaid on the measured differential 
cross section in Fig. 21 Also shown are results from 
an a ppr oximate next-to-NLO (NNLO) pQCD calcula- 
tion [2^ computed using MSTW2008 NLO PDFs 



271 ] calculations obtained 

nit = 



and same scales choices as the NLO result, and from 
the MC@NLO im (using CTEQ61 PDFs), ALPGEN, and 
PYTHIA event generators. The QCD scale uncertainty 
was evaluated for the NLO pQCD calculation [2^ 27 1 



FIG. 3: Comparison between the background-subtracted re- 
constructed top-quark pr spectrum and the one corrected for 
the effects of finite experimental resolution (two entries per 
event). Inner and outer error bars represent the statistical 
and total (statistical and systematic added in quadrature) 
uncertainties, respectively. 



by varying /i = mt = 170GeV/c^ by factors of 2 and 
1/2, and the PDF uncertainty by the approximate NNLO 
code [2§]. The total uncertainty is < 4% with only a 
small (< 1%) shape variation. A comparison of the ratio 
of da/dpT relative to a NLO pQCD calculation is shown 
in Fig. [5] The NLO pQCD calculations agree with the 
measured cross section, however, results from ALPGEN 
(pythia) have a normalization shift of about 45% (30%) 
with respect to data. A shape comparison of the ratio of 
(l/cr) da/dpT relative to NLO pQCD is shown in Fig.|6l 
All of the calculations reproduce the observed shape. The 
and corresponding probabilities [32| for the com- 
parisons in Figs. [5] and [6] of predictions to data are given 
in Table HVl 



In conclusion, we have presented a 1 fb~ measure- 
ment of the differential cross section of the top-quark 
Pt for tt production in the £-|-jets channel using pp col- 
lisions at = 1.96 TeV. Results from NLO and NNLO 
pQCD calculations and from the MC@nlo event genera- 
tor agree with the normalization and shape of the mea- 
sured cross section. Results from ALPGEN-f pythia and 
PYTHIA describe the shape of the data distribution, but 
not its normalization. We thank the staffs at Fermilab 
and collaborating institutions, and acknowledge support 
from the DOE and NSF (USA); CEA and CNRS/IN2P3 
(France); FASI, Rosatom and RFBR (Russia); CNPq, 
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FIG. 4: Inclusive da/dpr for production (two entries per 
event) in data (points) compared with expectations from NLO 
pQCD (solid lines), from an approximate NNLO pQCD cal- 
culation, and for several event generators (dashed and dot- 
dashed lines). The gray band encompasses uncertainties on 
the pQCD scale and parton distribution functions. Inner and 
outer error bars represent the statistical and total (statistical 
and systematic added in quadrature) uncertainties, respec- 
tively. 



FIG. 5: Ratio of dcr/dpT for top quarks in ti production (two 
entries per event) to the expectation from NLO pQCD. The 
gray band encompasses uncertainties on the scale of pQCD 
and parton distribution functions. Also shown are ratios rel- 
ative to NLO pQCD for an approximate NNLO pQCD calcu- 
lation and of predictions for several event generators. Inner 
and outer error bars represent statistical and total (statistical 
and systematic added in quadrature) uncertainties, respec- 
tively. 
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TABLE III: Inclusive differential cross section da/dpT for tt 
production at ^ = 1.96 TeV and m* = 170 GeV/c^ There 
are two entries per event, with the total normalized to the 
ti production cross section. In addition to total systematic 
uncertainties on the shape in pr in each bin, there is a pr- 
independent systematic uncertainty of 9.6% that is not in- 
cluded in the table. 
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TABLE IV: The x^/NDF and probability for comparisons 



PT 


(pt) 


Gross 
Section 


Stat. 
Unc. 


Shape Sys. 
Unc. 


NLO 
pQGD 


between the measured data and predictions using correlated 
(uncorrelated) uncertainties for the absolute (shape) compar- 


(GeV/c) 


(GeV/c) 


(fb/GeV) (fb/GeV) 


(fb/GeV) 


(fb/GeV) 


ison. 








0-45 
45-90 


29 
68 


70 
130 


11 
20 


5 
10 


59.6 
116 


Prediction 


Absolute 
x7NDF prob. 


Shape 
x7NDF prob. 


90-140 


113 


89 


13 


6 


83.8 


NLO pQGD 


0.695 
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0.315 0.904 


140-200 


165 


37 


6 


3 


35.6 


Approx. NNLO pQGD 0.521 
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0.497 0.779 


200-300 
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7.72 


mcQnlo 


1.22 


0.295 


0.777 0.566 
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329 


1.1 


0.3 


0.1 


0.814 


PYTHIA 


2.61 


0.0157 


0.352 0.881 


o-tt (pb) 




8.31 


1.28 




7.54 


ALPGEN 


5.04 


3.54x10"^ 


0.204 0.961 
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FIG. 6: Ratio of (1/cr) da/dpr for top quarks in tt production 
(two entries per event) to the expectation from NLO pQCD. 
The gray band encompasses uncertainties on the scale of 
pQCD and parton distribution functions. Also shown are ra- 
tios relative to NLO pQCD for an approximate NNLO pQCD 
calculation and of predictions for several event generators. 
Inner and outer error bars represent statistical and total (sta- 
tistical and systematic added in quadrature) uncertainties, 
respectively. 
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